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A better large N expansion for chiral Yukawa models.* 

George Bathas a and Herbert Neuberger a ' 

a Department of Physics and Astronomy 
Rutgers University, Piscataway, NJ 08855, U.S.A. 

We consider the most general renormalizable chiral Yukawa model with SU(3) co \ OT replaced by SU(N C ), 5C/(2)l 
replaced by SU(N W ) and U(1)y replaced by [/(l)^™ -1 in the limit N c — » oo, N w — » oo with the ratio p — » / jf- 7^ 

0, 00 held fixed. Since for N w > 3 only one renormalizable Yukawa coupling per family exists and there is no 
mixing between families the limit is appropriate for the description of the effects of a heavy top quark when 
all the other fermions are taken to be massless. The large N = V 'N C N W expansion is expected to be no worse 
quantitatively in this model that in the purely scalar case and the N = 00 limit is soluble even when the model 
is regularized non-perturbatively. A rough estimate of the triviality bound on the Yukawa coupling is equivalent 
to rat < 1 TeV. 
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We wish to study non-perturbatively a model 
that is reasonably close to what the minimal 
standard model would look like if the top quark 
Yukawa coupling and the Higgs selfcoupling were 
both large while all other couplings were kept at 
their known small values and could therefore be 
neglected in a first approximation. 

We propose to do this using an -^ expansion on 
a model defined non-perturbatively with a Pauli 
Villars type of regularizationpl]. The action is: 



S = 



5o 



Ao 



,2 ii* ±i 



ri* ±i\2 



4/V u 



T^y^x^ + x+W]; so >o (1) 

V N 



where 



N = VN W N C , 0+=a^, 



^1,2,3 = <7l,2,3j CT 4 = -(74 = —i 

i = 1, • • • , N w a — 1, • • • , N c 



(2) 



The 0X2,3 are Pauli matrices. Ultraviolet diver- 
gences can be eliminated by just regularizing the 
X and 4> propagators because there are no loops 
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consisting of only i\> propagators. This is achieved 
with h x {p 2 ) = h^p 2 ) = 1 + (p 2 /A 2 )" with n > 3. 
i\) represents the top-bottom doublet, x the 
right handed top and <f) the Higgs doublet. The 
index i generalizes the weak isospin index and 
the index a labels colors. ip+ is a Euclidean left 
handed Weyl spinor, X- a right handed one and 
<j> a complex scalar. ^_ and x+ are independent 
fields of handedness opposite to their unbarred 
counterparts as usual in the Euclidean functional 
integral formulation. 

The large N limit is taken with p = \ -jff- 
held fixed and is dominated by planar diagrams. 
Renormalizability allows us to get away with not 
having to introduce four Fermion interactions and 
with them an insoluble planar diagram problem, 
at no cost to generality]^ Thus the model is man- 
ageable being similar diagrammatically to two di- 
mensional QCD with matter at large JV C . The 
large N limit can be found by looking at the struc- 
ture of planar Feynman graphs: essentially one 
sums over all "cactuses of bubbles" and all "rain- 
bows" . For future work it is somewhat preferable 
to use functional integral manipulations to the 
same end. 

As an example of what can be done in this 
scheme we present a very simple calculation: 
We compute a /3-function associated with the 



1 Except for possible cutoff effects. 



Yukawa coupling constant defined in some sim- 
ple manner and ignoring cutoff effects. This can 
be done with relative ease because even at infinite 
order in the couplings the flow of this Yukawa cou- 
pling is unaffected by Ao, similarly to the finite N 
form at one loop order. From this calculation a 
rough estimate for the "triviality" bound on the 
top mass can be obtained.^ 

To find (3 functions it suffices to solve the model 
at criticality when all masses vanish. We start in 
the symmetric phase and pick a mass independent 
renormalization scheme to make the zero mass 
limit smooth. 

We introduce bilocal fields, 



S a ,p( x >y) = jfX a + a ( x )x a -p{y) 



(3) 



where a, (3 are space-time spinor indices, into the 
functional integral by writing representations for 
the appropriate ^-functions with the help of two 
auxiliary bilocal fields \ a> p(x,y) and p(x,y). We 
can integrate out now all the original fields, in- 
cluding the zero mode of the scalar field because 
we assumed that we are in the symmetric phase. 
This makes the TV and p dependence explicit and, 
when N — > oo, the integral over S,X,K and (i is 
dominated by a saddle point which has to satisfy 
the following equations: 

K(x,y) = {[-h$(-d 2 )d 2 + ml + A*]" 1 },,,* 

S a ,p{x, y) = {[h x (-d 2 )0- + A] } y 0, xa 

X a ,p(x,y) = -g 2 pV2{0^ 1 } xa>y0 K(x, y) 

Q 2 
ftfay) = ■z-7=z2{0+ 1 }x a , v 0S a ^(x,y) 



pV2 



a,0 



+^-K(x,x)5 Xty (4) 



We go to Fourier space and use translational and 
rotational invariance to write 



Aq,/3(p) = ip^o-ftXip 2 



jS(p) = p(p 2 ) (5) 



Let us introduce Z x and Z, 



wave function renor- 
malization constants, and new functions f(p) = 



2 In practice this bound might be somewhat uninteresting 
because it is likely to be higher than the vacuum stability 
bound for reasonable Higgs masses. 



Z x [\{p 2 ) + h x {p 2 )) and g{p 2 ) = Z^h^p 2 ) + 
Wq + p{p 2 )]. Since at the planar level there are 
no vertex corrections to the Yukawa coupling]^] it 
makes sense to define a renormalized Yukawa cou- 

2 

pling by g% = Z^Z x g 2 , (a = I |f T ). We also set 
Ag = AoZ? but, unlike a, do not expect Aq to stay 
finite when A — > oo. Z x and Zq are defined in 

terms of the bare couplings by f(p 2 ) = ^ 2 = 1< 
After performing the angular part of the mo- 
mentum integrals one gets when the products 
in the saddle point equations for x,y dependent 
quantities are converted into convolutions in mo- 
mentum space one finds: 



f(p 2 ) = Z x h x (p 2 ) + V2ap 
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The massless case is obtained by adjusting m 2 , so 



that 



= and then it makes sense to define 



h(p 2 ) — 9 <z ' leading to a set of equations of a 
more symmetric appearance with an interesting 
bosonic-fermionic symmetry at p — 1 . 

Rescaling the functions and the momentum 
variable by 



/(V) = 

we obtain: 



2a<j)(u) h(up ) = \faip{u) 



(7) 
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^(«) = -7= + -tMM u M 2 ) - h<t,(p 2 )} 



1 f 1 dv f 1 (1 - x)dx 
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(«) = irr + ;rHM«/* 2 ) - h x (p 2 )] 

2y/a 2y/a 

_ f 1 dv f 1 (1 - x)dx 



u v J xp{xv) 



(8) 



'This is also true at one loop order at finite N. 



It is clear that when n is large enough h and 
/ are dominated by the /i^. x functions in the ul- 
traviolet and therefore the large cutoff behavior 
of the wave function renormalization constants is 
such that in the infinite cutoff limit, at u and p 2 
fixed, Z,), and Z x simply disappear from the above 
equation, leaving us with 

f 1 dv f 1 (1 - x)dx 

iol P Ju v 

1 f 1 dv . , , 

P — i, \ ( 9 ) 



ip(u) 



o Hxv) 
1 (1 - x)dx 



2 V" r Ju v Jo l/j{xv) 

On physical grounds it is obvious that we wish 
that ip and <f> be positive. The equations above 
show that if this is true in some interval (0, m*) 
then the functions will be monotonically decreas- 
ing there. However, the equation do not admit an 
asymptotic behavior as u — * oo with both func- 
tions approaching non-negative limits so the pos- 
itivity requirement must get violated somewhere 
in the ultraviolet Jj Starting at the scale where 
the violation first occurs cutoff effects cannot be 
neglected any more. This "unphysical" scale is 
the usual Landau pole, this time appearing in a 
nonperturbative approximation. 

In the infrared there is no Landau pole problem 
and cutoff effects are indeed negligible. To inves- 
tigate the behavior there we set x — — log(u) and 
ijj(x) — V'(m), 4>{x) = 4>(u) and derive: 



V 



p<P 

d f d 



dx V dx 



V4> = Z 
d 

dx 



(10) 



These equations allow us to extract the infrared 
behavior of the solution of the integral equa- 
tions. We ended up with only differential equa- 
tions (rather than integral) reflecting that accord- 
ing to the Renormalization Group one needs only 
very limited information at a given scale in order 
to derive the behavior at a scale close by; there- 
fore a differential equation must show up even- 
tually, its order less the number of asymptotic 
conditions determining the amount of informa- 
tion at a given scale that is needed to determine 

4 A zero of i/5 or <f> at a positive u corresponds to a pole 
in a two-point function at an Euclidean momentum; the 
"particle" associated with this pole would be tachyonic. 



the behavior at the next scale. Our choice of reg- 
ularization was made so that even at finite cutoff 
<p and ip obey a set of purely differential equations 
of a structure similar to |l0]. 

Since the fixed point governing the infrared be- 
havior is the free field fixed point the form of the 
solutions in the infrared simply embodies the two 
anomalous dimensions associated with the x an d 
<p fields when expanded in a around a = 0. These 
anomalous dimensions as a function of the cou- 
pling also determine the /3 function and from the 
asymptotic series of ip{x) and 4>{x) at x — v oo 
we can get the contributions to the (3(a) function 
ordered in the number of loops. 

The differential equations lead to: 



1 i &n i logx i ' l i i 

ip — yjox [1 + c h c - + ■■ ■] 

x x 

$ ~ o z a [l + rf^ + d'~ + ■ ■ ■} 

X X 

where 

,-2 P 2 , -2 3 P 

a = op = rr, d = cp = 



(11) 



d' 



1 + p 2 
3(1 -P 2 ) 
2(1 + p 2 ) 2 



^000 



(1 + p 2 ) 3 
P + P- 1 



(12) 



Only the product ipofio and a particular linear 
combination of d! and c' get determined because 
the equations are invariant under a field rescal- 
ing ip —> A~ lr i/j, <p ~^ A(p reflecting a change in 
the finite parts of the wave function renormaliza- 
tion constants and a shift x — > x + xa reflecting a 
change in p? . The field-rescaling invariance dis- 
appears in the product r(u) = y/2ip(u)<f>(u) and 
this combination is indeed special because it sat- 
isfies a Renormalization Group equation in which 
only the /3-function appears (the sum of the two 
anomalous dimensions rather than each individ- 
ual one): 






r{e 2 \a) =0 



(13) 



This equation is exact in our limit. From it we 
derive another exact relation between r and j3: 



/3(a(V)) = -2[ a (V)] 2 ^T 



(14) 
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Figure 1 . The nonperturbatively determined /?— 
function for the Yukawa coupling at N c , N w — > oo 
with M*- = | . At tree level the top quark mass is 
given by m t = (47r 2 a/3)i 246 [GeV] : /3(a) = ^ 
where i is the logarithm of the energy scale. 



Using now the asymptotic expansion of the ex- 
act solution in the infrared we get: 



(3(a) = a 2 [V2(p + p- 1 ) - 3a + 0{a 2 )] 



(15) 
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^s which is close to the finite N result, 
The difference comes from the wave 



The one loop result when expressed in a con- 
ventional variable, y(t) (in the standard model 
the top mass at tree level is m t — j/(~ 
0) 246 [GeV]), is, with the true value p — 

dy 2 _ 
dt 
dy 2 _ 9y 

dt 8;r 

function renormalization of the left handed ip field 
which is suppressed at leading order in i . For ar- 
bitrary N c and N w we would have obtained at one 
loop order -J^- = - — c+ 87r ™ + and we see that 
the relative magnitude of the 1/N correction is 
much smaller than in the scalar case||. 

To get the full /3-function at infinite N we solve 
for ip and <fi numerically by iterating an integral 
form of the equations @ but with the normaliza- 



tion conditions at p? incorporated; this avoids the 
appearance of the double integrals in [8| and |9| . 

The result for /3 is shown in the figure and 
one again sees explicit evidence for the nonper- 
turbative existence of a Landau pole. We com- 
pare the Landau pole energy obtained from the 
exact result to the one loop estimate. For exam- 
ple, with the high value of a(p 2 ) = 1.365 chosen 
in the figure we obtain A 2 jandau = 2/x 2 while at 
one loop A 2 andau = 1.565/U 2 . A top quark cor- 
responding to such a strong coupling would have 
a mass of m t = ^^ 246 [GeV] ~ 1043 GeV. 
We see that "triviality bounds" on the top will 
likely come out close to the perturbative unitarity 
bounds of Chanowitz et. al.|| in contrast to Ein- 
horn and GoldbergB who obtained 5600 GeV. 

From our non-perturbative result in the figure 
one sees that the one and two loop errors in (3 are 
about comparable for a w 0.4 (roughly of order 
20 percent) and that using the one loop results all 
the way out to infinity is not very harmful, but the 
two loop result would be completely misleading 
if extrapolated beyond its region of validity, to 
a«l. A top mass of 200 GeV would correspond 
to a « 0.05, well within the perturbative domain. 

For the future we plan the following steps: 

• Compute the /3-function associated with A at 
infinite cutoff in the massless case. 

• Reintroduce the finite cutoff in the above 
computation to see the set in and magnitude of 
cutoff scaling violations. 

• Redo both the infinite cutoff and the finite 
cutoff computation away from the critical line, i.e. 
when thresholds are present. This has to be done 
in both the broken and the symmetric phases. 

• Compute the effective potential at infinite 
cutoff limit first and then at finite cutoff and ana- 
lyze the issue of "vacuum instability" in this non- 
perturbative setting. 

• Study bound state problems. 
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